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Abstract 



Using e + e~ — > and e + e~ — > qq(7) events radiative to the Z pole, 

DELPHI has determined the centre-of-mass energy, a/s, using energy and mo- 
mentum constraint methods. The results are expressed as deviations from the 
nominal Lep centre-of-mass energy, measured using other techniques. The re- 
sults are found to be compatible with the Lep Energy Working Group estimates 
for a combination of the 1997 to 2000 data sets. 
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1 Introduction 

One of the main goals of Lep2 is the precise measurement of the mass of the W-boson, 
Mw- The aim, with the integrated luminosity collected, is a total experimental uncer- 
tainty of about 40 MeV/c 2 . One important contribution to the systematic uncertainty 
on M w comes from the accuracy with which the centre-of-mass energy, E CM , can be 
measured. To match other statistical and systematic errors, the error on E CM at Lep2 
should be below 30 MeV, equivalent to an error on the beam energy of less than 15 MeV. 
At LepI, very accurate determinations of Ecm were made which relied on precise beam 
energy (Eb ca m) measurements provided by the resonant depolarisation technique [1]. At 
Lep2 energies resonant depolarisation is no longer possible. Therefore a method using 
a magnetic extrapolation is employed, in which resonant depolarisation is used at beam 
energies of 44-61 GeV to calibrate NMR probes in selected dipole magnets. These probes 
are then used to estimate the beam energy at Lep2 physics energies by extrapolating 
from the lower calibration energies [2]. 

The uncertainty on the extrapolation procedure dominates the energy determination. 
Several methods have been developed by the Lep Energy Working Group to constrain this 
uncertainty [2]. These include the measurement of the beam deflection in a spectrometer, 
an analysis based on the evolution of the synchrotron tune with the RF voltage, and a 
comparison of the bending field as measured by the NMRs with that sensed in the ring 
as a whole by the flux loop. The assigned error on E CM from combining the results 
of these analyses is 20-25 MeV for the majority of Lep2 operation, with a larger error 
of about 40 MeV in 2000, where in addition to the magnetic extrapolation itself, other 
uncertainties are significant [2]. 

The collision energy can also be measured in both hadronic and leptonic channels 
using the "radiative return" events. In these events the effective e + e~ centre-of-mass 
energy is reduced to that of an on-shell Z-boson by initial-state radiation (ISR). The 
Aleph collaboration performed a measurement based upon the latter method using just 
the 183 GeV data [3]. Subsequently, measurements at this and higher energies have been 
carried out by the L3 [4] and Opal [5] collaborations. This note describes a method that 
uses e + e~ — > fi + fi~ (7) and e + e~ — > qq(7) radiative return events to measure the difference 
between the centre-of-mass energy determined from radiative events and the value from 
the Lep Energy Working Group AEqm- This method, particularly if results from the four 
LEP experiments can be combined, provides a cross-check against large systematic effects 
in the Lep Energy Working Group determination. The dimuon channel is the easiest in 
which to study the radiative return peak due to its simple two-prong topology. However, 
the statistical power is much smaller than that available from hadronic radiative returns. 

The analysis used the data of the four last years of Lep data taking with the collision 
energy ranging from 183 to 207 GeV. The total data sample was organized in 8 energy 
points: 183, 189, 192, 196, 200, 202, 205, 207 GeV. The measurement was performed at 
each energy point separately and at the end a combined result was produced. Details 
of the components of the Delphi detector and their performances can be found in [6,7]. 
Data from the period after 1st September 2000, when part of the TPC, the main tracking 
device of Delphi, had high voltage problems, were not included. 

Simulation of signal and all relevant backgrounds was available for each energy. The 
e + e~ — * /i + /i _ (7) and e + e~ — > qq(7) events were simulated with the KK 4.14 genera- 
tor [8]. Most of the four-fermion final states were produced with the program described 
in [9], based on the WPHACT 2.0 generator [10]. However, as discussed in [9], for the 
final states e + e~ff, collisions with either the e + or e~ emitted at less than 2° to the 
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beam, and with the invariant mass of the ff system less than 40 GeV/c 2 were generated 
with PYTHIA 6.143 [11]. Hadronisation was performed with PYTHIA 6.156. In addi- 
tion, ARIADNE 4.08 [12] was used for fragmentation studies. The generated events were 
passed through the DELPHI detector simulation program DELSIM [7], and then through 
the same reconstruction and analysis programs as the data. 

Section 2 describes the analysis of the e + e~ — > channel, and Section 3 that 

of the e + e~ — > qq(7) channel. Conclusions are given in Section 4. 



2 The e + e — > (7) radiative return events 

2.1 The determination of \fs i 

For a centre-of-mass energy, y/s, above the Z resonance the differential distribution 
of the invariant mass of the fj, + fx~ pair, y/s*, has two distinct peaks. One is a peak at 
y/s' Rj \fs and the other is at y/s' ps M z . The latter peak is formed by the radiation 
of one or more photons from the incoming electron or positron, reducing the effective 
invariant mass of the e + e~ system to that of an on-shell Z, where the cross-section for 
s-channel Z exchange is large. 

The lineshape at the Z has been studied in great detail at LepI. The measurements of 
M z and T z have been made to MeV precision [13]. Using the knowledge of the resonance, 
combined with that of QED initial-state radiation, one can hope to model the radiative- 
return peak to similar precision [14]. 



If a variable x = ys'/s is defined, a distribution of x will have some value, xz, at 
the position of the radiative-return peak. In a distribution of y/s*, the position of the 
radiative return peak is always at some value (\Js^ ~ M z ), independent of the initial 

centre-of-mass energy. Therefore, if ^Js^ is known, a measurement of xz can be used to 
calculate y/s from the expression: 



However, the measurement of xz must be independent of y/s for this to be a valid 
method of evaluating \/s. This is true for radiative fi + fi~ events, which are consistent 
with only one hard ISR photon being radiated collinear to the beams. 

For this simple case of planar events s' can be expressed in terms of only s and the 
energy of the radiated photon, E^. 

s' = s - 2E^. (1) 

Using 4-momentum conservation, E 1 can be written in terms of the polar angles (defined 
relative to the incoming electron beam) of the muons and ^/s alone: 

E |sin(^+ + flg-)]vg (2) 

7 ~" sin(V) + sin(^-) + |sin(V+^-)l' 

where 0^- and 0^+ are the polar angles of the momentum vectors of the muon and anti- 
union, respectively. Substituting for E 1 in equation (1), leads to an expression for x: 



= x = 



sin(6^+) + sin(6» M -) - 


sin(^+ + ^-) 


sin(0 At +) + sin(6» M -) + 


sin(V + V)l 



(3) 
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in terms of only 6> M + and 9^- , as desired. 

To investigate the possibility of a systematic discrepancy between y/s from the method 
and that from the standard Lep measurement, it is also useful to be able to combine 
data from many years. A quantity is defined that expresses the difference between the 
centre-of-mass energy by the method of using radiative-return events and the value from 
Lep, 

AE CM = v^Angles _ v^Lep = ~~ v^Lep- ( 4 ) 

If the two methods agree, this should give a distribution, centred around zero, that 
takes the form of the radiative-return peak. A fit can be performed to this distribution 
and a central value, giving the overall difference, extracted. 



2.2 Event selection 

Initially the standard selection for dimuon events, as used in the cross-section and 
asymmetry measurements, was applied. A detailed description of the selection can be 
found in [15]. To reject dimuon events which are not consistent with a single ISR photon 
down the beam pipe, further cuts were made: 

• it was required that no photons were seen within the electromagnetic calorimetry. 
For this study the calorimetry includes the low-angle luminosity detectors, STIC, as 
well as the barrel and forward regions, collectively called ECAL. Photons of energy 
greater than 500 MeV in the STIC, and greater than 2 GeV in the ECAL were 
considered. This cut removed events with a high transverse-energy photon; 

• the acoplanarity 1 of the two muons with respect to the beam axis was required to 
be less than 0.1 rad.; 

• the ratio of the absolute value of the missing longitudinal momentum (i.e. net 
momentum of the two muons along the electron beam direction) to the missing 
energy (i.e. the nominal Lep2 centre-of-mass energy minus the sum of the two 
muon energies) in the event had to be greater than 0.8. This selection cut is effective 
against multiple ISR photons. 

The latter two selections also reduce backgrounds from r + r~, four-fermion events and 
from cosmic rays. This leaves a background level of approximately 2% in the data set 
used in this study. 



2.3 Fitting the radiative return peak 

The selected events were fitted to the expected shape of the differential cross-section 
in AEcm using an unbinned maximum log-likelihood fit. The function, f(x), used is: 



where Af, x%, r z FIT , B\ and £>2 are parameters determined by fitting. This is an empirical 
function which fits well the shape of the simulated data. The package MINUIT [16] was 
used to maximise the likelihood function. 

Events generated by KK [8] and passed through the Delphi detector simulation pack- 
age were used to calibrate the fit. It was found that the central fitted value, calculated 

1 Thc acoplanarity is defined as the complement of the angle, in the plane transverse to the beam, between the two 
tracks. 
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from the parameter x%, did not correspond to 0, but was systematically shifted to lower 
values of AEcm, lower values of -y/i. A number of different samples at several centre- 
of-mass energies were used to evaluate the size of this bias. The bias was found to be 
independent of energy over the range of energies considered. Thus the bias from the 
combination of all the energies is used as the correction to the final fitted central value 
from data. 

The function (5) was fitted to DELSIM samples for KK events generated at centre- 
of-mass energies at or close to those of the data. The result of the fit to all simulated 
data is shown in Figure 1. The bias from this fit was found to be AEcm — —0.057 ± 
0.022 (stat.) GeV (x 2 /N.D.o.F. = 26.8/25). 

DELSIM 

•^10000 F~ 



5u 




-75 -10 -5 5 10 15 
AE CM (GeV) 

Figure 1: A fit to the AE CM distribution (see Section 2.1) of fully-simulated 
e + e~ — > /z + /i~(7) events using the KK generator [8], corresponding to an effective lumi- 
nosity of about 86 fb -1 . The points show the simulated data, and the shaded histogram 
is the fitted function. 



2.4 Results 

From the 1997-2000 data sets 600 events were selected and the AEcm distribution was 
fitted with equation (5), with the parameters £>i and £> 2 fixed at their values obtained 
from simulation. The fit results at each energy point are shown in Figure 2. Since, within 
the large statistical errors, there is no apparent dependence on energy, the 1997-2000 data 
sets were combined. The fit is shown in Figure 3. The result of the fit was: 

AE CM = +0.184 ± 0.150 (stat.) GeV (x 2 /N.D.o.F. = 8.6/12). 

The statistical error was verified using the "bootstrap" method [17]. In this method, 
random subsamples of simulated events of the same size as the real sample size are 
selected and for each of them the result is computed in the same way as for the real sample. 
This gives some distribution of the resulting value. The spread of this distribution was 
compatible with the statistical error assigned by MINUIT (i.e. the "pull value" was found 
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DELPHI 



0.5 - 



- 



-0.5 - 



-1 - 



180 185 190 195 200 205 210 

(GeV) 



Figure 2: Measured centre-of-mass energy shift at each energy point from 
e + e~ — > /i + /z~(7) events. The error bars show the size of the statistical error. The line 
represents the fit to all energies combined. 



DELPHI 




CM 



10 15 
(GeV) 



Figure 3: AEqm value for data (points) corresponding to the combined sample of 
e + e~ — > /i + /i~(7) events collected from 1997-2000. The fitted function is shown by the 
light grey histogram. The dark grey histogram shows the expected background. The 
integrated luminosity is 594 pb _1 . 



to be 1.01). However the average statistical uncertainty on the centre-of-mass energy of 
the samples was found to be 0.172 GeV, somewhat larger than that found with the real 
sample. 
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2.5 Systematics 

The systematic uncertainties are presented below. They include contributions from 
the modelling of the ISR spectrum as well as effects expected from reconstruction and 
from the Lep beam parameters. 

2.5.1 Systematics from the bias 

The bias estimation has a systematic uncertainty on AE CM of ±22 MeV, due to the 
limited statistics of fully-simulated events. 

2.5.2 Systematics from radiative photonic correction modelling 

The effects of the QED modelling were studied at generator level using the KK [8] 
program. The KK generator calculates photon emission from initial beams and outgoing 
fermions to second order, and includes ISR®FSR interference also to this order. 

To make some estimate of the magnitude of the effect of imperfections in the radiative 
photon modelling on the bias, samples of 500,000 KK events were generated at each 
of the centre-of mass energies used and fitted with the function given in equation (5). 
Following the advice given by the authors of KK [8] , an estimate of the uncertainty of the 
ISR modelling was obtained by weighting the events generated with the default CEEX2 
(0(a 2 )) approximation to correspond to the CEEX1 (0(a)) approximation, and fitting 
equation (5) to the weighted events. Half the difference in the resulting bias was taken 
as the systematic uncertainty due to this effect. It amounted to 2.5 ±0.7 MeV on AEcm- 

2.5.3 Mismeasurement of the DELPHI aspect ratio 

The aspect ratio is defined as the ratio of the length to the width of the detector. It is 
limited to the precision to which the position and dimensions of the Vertex Detector (VD) 
can be measured. The effect of a mismeasurement of the aspect ratio is to introduce a 
bias on the measurement of the polar angle, 6. A positive (negative) bias on tan# would 
systematically increase (decrease) the opening angle of the two muons and lead to an 
underestimate (overestimate) of y/s. 

The correspondence of hits in overlapping silicon modules is sensitive to a misalignment 
of the Vertex Detector. In fact the study of these overlaps constitutes an essential part of 
the procedure for the alignment of the Vertex Detector. The alignment was carried out 
independently for each year of operation. A careful analysis of the data of all years shows 
that a residual bias of 6 x 10~ 5 can be excluded at the 95% confidence level. However a 
study of the behaviour of the Delphi silicon modules [18] shows that the centre of charge 
of the holes (and electrons) differs from the mid-plane of the detectors by 10 — 20 /im. In 
the alignment procedure this difference would be automatically accounted for. However 
much less is known about the displacement of the electrons (used in the measurement 
of z, along the beam axis, and hence 6) since the necessary overlaps are much fewer. 
Taking ±10 /im as an effective uncertainty on the location of the sensitive region of the 
detectors, and allowing for the fact that the polar angles of some tracks are determined 
from signals in the silicon detectors themselves (without using the interaction point as a 
constraint) it is concluded that a reasonable estimate of the aspect ratio uncertainty is 
3 x 10~ 4 . Such a bias would correspond to a shift of AE C m of about ±29 MeV. 
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2.5.4 Systematics from LEP beam parameters 

In general the colliding e + and e~ are not collinear at the interaction point, and their 
energies are not equal. This leads to a boost of the centre-of-mass system relative to 
Delphi and in an individual event will lead to a change in the fitted centre-of-mass 
energy. However, averaged over all event orientations, the effects are much reduced. 
Similarly, inequality of the beam energies or a tilt of the beams with respect to the 
Delphi coordinate system will lead to errors in the fitted energies. Studies at generator 
level with beam acollinearities and tilts up to 500 yurad and beam energy differences up to 
100 MeV, have shown that the overall effects would be approximately 1 MeV. However 
the spread in energies of the electron and positron beams, which can give an rms spread in 
the centre-of-mass energy of ±250 MeV at high energy, results in a correction of 6 MeV. 

2.5.5 Muon polar angle and momentum reconstruction 

Studies at generator level showed that random errors in the muon polar angle due to 
reconstruction inaccuracies in DELPHI lead to systematic shifts of the fitted energy. In 
principle, if correctly estimated, these are taken into account in the simulation. The selec- 
tion on energy and momentum imbalance also makes use of the measured momenta of the 
muons, so in principle is affected by bias or additional smearing in this quantity. However 
it was estimated that the overall effect of these imperfections led to an uncertainty in the 
centre-of-mass energy of less than 2 MeV. 

2.5.6 Backgrounds 

Backgrounds from four-fermion and r-pair production have been included in the sim- 
ulated samples from which the bias was calculated. Inclusion of the 2% estimated back- 
grounds changed the bias by 0.4 MeV on AEcm 5 with a statistical uncertainty of ±4 MeV, 
due to the size of the samples of simulated background events. The latter value was taken 
as the systematic uncertainty due to the background. The energies of the simulated back- 
grounds were not varied during the fit. It was estimated that this introduces a negligible 
additional systematic uncertainty. 

2.5.7 Knowledge of the Z mass 

As can be seen from equation (4), there will be an uncertainty in the estimate of 
AEcm, due to the uncertainty of the Z mass as measured at LepI [13]. The resulting 
uncertainty from this source amounts to ±4 MeV. 

2.5.8 Summary of systematic uncertainties 

Table 1 gives a breakdown of the systematics on Eqm- The total systematic uncertainty 
on AE CM is 38 MeV. Thus the dimuon events give a value of: 

AE CM = +0.241 ± 0.150 (stat.) ± 0.038 (syst.) GeV. 
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Source 


Hvrror on AUvcm (,Mev ) 


Bias 


22 


0"Fj"D m nnpl li n q* 


3 


Aspect ratio 


29 


Lep parameters 


6 


Angular resolution 


2 


Backgrounds 


4 


Z mass 


4 


Total systematic error 


38 


Statistical error 


150 



Table 1: The sources and associated values of uncertainties on AEcm, measured from 
radiative dimuon events. 

3 The e + e~ — > qq(7) radiative return events 

3.1 Track and event selection 

Events of the type e + e~ — > qq(7) give rise to a large hadronic multiplicity. The anal- 
ysis was based mainly on the charged particle tracks reconstructed in each event. Track 
selection criteria were applied in order to reject tracks originating from secondary inter- 
actions as well as poorly measured tracks. The charged particles were required to have 
momentum larger than 0.4GeV/c, but not exceeding 1.5 times the beam energy, momen- 
tum resolution better than 100%, and transverse and longitudinal impact parameters less 
than 4 cm and 10 cm, respectively. A more detailed description of the selection can be 
found in [15]. 

All hadronic events were required to have at least 7 charged particles and a transverse 
energy greater than 20% of the collision energy. The transverse energy was calculated on 
the basis of the transverse momentum of each selected track and on the electromagnetic 
showers reconstructed in the HPC and the FEMC with a shower energy above 500 MeV. 
No other cuts were used, in particular for suppression of the WW and ZZ backgrounds, in 
order to avoid a possible bias to the results. The selection efficiency for hadronic events 
was slightly less than 80%. 

3.2 Fitting of the radiative return peak 

In the s' determination a two-jet configuration was forced. If no ISR photon candidate 
was observed then the reduced energy was computed from the angles of the jets assuming 
a photon escaped into the beam pipe. If an isolated energetic photon was observed then 
the value of s' was derived from the combined jet and photon directions, assuming that 
additional energy could be radiated inside the beam pipe. If no solution was found, the 
reduced energy was derived from the measured energy of the isolated photon alone. 

Figure 4 compares \fs' distributions between data and simulation [8] around the ra- 
diative return peak. In the analysis the interval of 75-105 GeV of the a/s 7 distribution 
was used with bins of width 1 GeV. 
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Figure 4: Reduced energy distribution of e + e~ — > qq(7) events. Points represent the 
data, light grey area the signal simulation, dark grey area the background simulation. 
The integrated luminosity is approximately 600 pb" 1 . The distributions combine all 
energy points 183-207 GeV. 
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3.2.1 Method of Monte Carlo reweighting 

The method of this analysis was to compare the Z-boson mass, M z , measured from 
data, to its nominal value and to find the centre-of-mass energy shift, AE CM , from their 
difference, AM Z , using: 

AM Z 

AE CM = — — — • E CM . 

1V1 Z 

The approach used for the M z determination is based on varying the M z value in 
simulation and finding the value which provides the best agreement of the y/s* distribution 
between data and simulation. M z was varied by reweighting the simulation using a Breit- 
Wigner-like function f(s',M z ): 

wetght(s',M z new ) /( ' 



f(s',M z ) 



f(s',M z nom y 
1 + P- (Vs'-Mz)) -s' 



(s'-M z 2 )2 + ( s '.r z /M z )2' 

where P is a skewness parameter (set to 0.01 GeV -1 ), yfs 1 is the generated reduced energy 
in the given event, M z nom the nominal value of M z and M z new the new value of M z . This 
particular form of the function /(«', M z ) was chosen because it describes well the shape 
of the Z mass peak of the generated y/s 7 distribution, and this is used to set the value of 
P. 



3.2.2 Fitting simulation to data 

The determination of M z was performed by fitting simulation to data using a log- 
likelihood method. The log-likelihood of the real distribution, given the simulated y/s' 
distribution, was computed as: 

n 

C = 2j2(fi-yikyfi), (6) 

i=i 

where n is the number of bins in the distribution histogram, and fi are the content of 
the i-th bin in the real and simulated distributions, respectively. 

For a cross-check a chi-square method was used as well. The chi-square of the difference 
was computed as: 

fi Hi \ 



fy% 

When the simulation has been reweighted to a new M z value the content of the simu- 
lated distribution fi...f n was recomputed accordingly. By a detailed scan, the dependence 
of x 2 and C on M z was found to be quadratic over the region of interest. It can be written 
as: 

2 

i 

Three points are needed to find the minimum of the parabola (8). The most convenient 
choice of points is M z — 0.5 GeV, M z and M z + 0.5 GeV. After reweighting the simulation 
to M z — 0.5 GeV and M z + 0.5 GeV and computing C or x 2 at all three points it is 
straightforward to compute M z mm as well as a, which gives the statistical error on M z mm . 




, 2(Mz) = w + ( M --y'"> 2 . ( 8 ) 
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3.3 Corrections in the forward region 

One of the main problems to be dealt with in the analysis was a noticeable excess of 
tracks at low polar angles (forward tracks) in data as compared to simulation. The effect 
is shown in Figure 5a. The most likely cause of this effect is an underestimation in the 
simulation of the track reconstruction efficiency for low-momentum particles at low polar 
angles. 

A tuning was applied to the data in order to correct this discrepancy. It consisted in 
removing some fraction of forward tracks on a random basis. A track could be removed 
from an event with some probability V computed from a ratio between the number of 
simulated N sim and real N rea i tracks: 

y ly sim 

N r eal 

In computing the correction it was taken into account that the size of the discrepancy 
depends on the track momentum and polar angle. The momentum range, as well as polar 
angle range, were split into a certain number of bins. The ratio N sim /N rea[ was computed 
in each bin and the correction was applied accordingly. 

For computation of the correction ratios N sim /N rea i, the data collected in calibration 
runs at the Z peak were used. Such data are available for every year of Lep2 data taking. 
They provide much more statistical power for computation of the tuning coefficients than 
the high energy data, and the corrections are uncorrelated with the Lep2 beam energy. 

Figure 5b shows the polar angle distribution after the correction is applied. Some 
residual difference still remains because the size of the discrepancy in high energy and in 
on-peak data was slightly different. 

This correction changed significantly the results because radiative events are boosted 
forward and therefore are sensitive to changes in forward tracking. The measured average 
shift in the centre-of-mass energy changed as a result of the correction from +206 MeV 
to -116 MeV. 

It should be noted that this effect does not influence the muon-pair result. The defect 
in the simulation would show itself as a difference in efficiency between simulation and 
data. However as part of the dimuon analysis, the detection efficency in simulation is 
tuned to that measured in data. Furthermore the effect is concentrated at low track 
momenta. Hence the results presented in Section 2 are not affected. 

3.4 Results 

The results obtained from fitting simulation to data by the log-likelihood method at 
each energy point are shown in Figure 6. The average shift of the measured centre-of-mass 
energy with respect to the Lep value is found to be —116 ± 106 (stat.) MeV. 

In the analysis the statistical error on M z mm is computed as a from equation (8). In 
order to verify that it gives a correct estimate of the statistical error, two other methods 
of the statistical error determination were used. 

The first one is the so-called "Jackknife" method [19]. It describes the measurement 
of the change in the result if one single event is taken out of the sample. If one event is 
removed out of the bin i, the change of the result R is A; = R — Ri. The variance a R 2 
of the result is then equal to the sum of A; 2 over all bins times the number of entries y i: 
or = J2i Aj 2 • Thus gr should give the same error on M z as before. 

The statistical error was also checked using the "bootstrap" method [17] described in 
Section 2.4. 
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Figure 5: Distribution of track polar angle of radiative return e + e~ — > qq(7) events (a) 
before and (b) after corrections in the forward region. Points represent data, light grey 
area the signal simulation, dark grey area the background simulation. The distributions 
combine all energy points 183-207 GeV. 



Statistical errors computed with the methods described above were found to be con- 
sistent within a few percent with the error computed from equation (8). 

The contribution of the statistics of the Monte Carlo sample amounts to an uncertainty 
of ±14 MeV on AE CM . 

3.5 Systematics 

3.5.1 Uncertainty on corrections in the forward region 

The stability of the Ecm determination was studied using alternative methods for the 
determination of \/F (using also neutrals or a kinematic fit). A correction procedure 
based on the high energy data rather than calibration data at the Z peak was also tried, 
although this could not be used for a result since the correction is not independent of the 
assumed centre-of-mass energy. A reweighting of tracks in data, rather than their removal, 
was also used. Considering the spread of the results, and also the size of remaining data- 
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Figure 6: Measured centre-of-mass energy shift at each energy point from e + e — ■> qq(7) 
events. The error bars show the size of the statistical error. The line represents the 
average (x 2 / N - D -°- F - = 9.4/7). 



simulation discrepancies (for example those visible in Figure 5b) it was estimated that 
a systematic uncertainty of ±120 MeV on Ecm should be assigned, due to the forward 
track correction and other imperfections in the simulation. 

3.5.2 ISR modelling uncertainty 

To study the uncertainty related to ISR modelling in simulation, two e + e~ — > qq(7) 
Monte Carlo samples, made with different generators, were compared. One sample was 
generated with KK and another one with PYTHIA 6.153. 

The comparison revealed no statistically significant effect. When one simulation was 
fitted with the other, the deviation of M z from the nominal value was found to be con- 
sistent with zero. The statistical error of this measurement was approximately 20 MeV 
on E C m- 

In addition, as for the dimuon events, studies were made at generator level using the 
KK generator. Samples of 100,000 e + e~ — > qq(7) events were generated at each of the 
energies used. Half of the difference between the CEEX2 and CEEX1 approximations 
amounted to 0.2 ± 0.4 MeV on AE CM . However in the DELPHI simulation, events are 
generated with KK without ISR0FSR interference. The hadronisation and FSR were 
generated using PYTHIA. The effect of removing ISR&FSR in KK-generated events 
amounted to a change in AEqm of 4.4 ± 5.4 MeV. 

The uncertainty connected with the ISR modelling was therefore taken as ±6 MeV on 
AE CM - 

3.5.3 Fragmentation uncertainty 

A similar approach was used for estimating the fragmentation uncertainty. A sample 
generated with KK and hadronised with PYTHIA was compared to a sample generated 
with KK and hadronised with ARIADNE [12]. Both PYTHIA and ARIADNE had been 
consistently tuned to the Delphi data. Fitting one simulation to the other gave a few 
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MeV shift in M z , with the statistical error on the difference of about lOMeV. Therefore 
the fragmentation uncertainty was estimated to be ±20 MeV on E CM - 

3.5.4 Four-fermion background uncertainty 

The majority of the four-fermion background is concentrated at high values of s'. 
Around the Z mass peak the background amounts to 8-10% (depending on energy point) 
of the signal and it has an almost flat distribution with s' (see Figure 4). Thus, a bias in 
the background will affect the results only weakly. 

The uncertainty due to four-fermion background was estimated by scaling the back- 
ground cross-sections by ±5%. It changed the result by ±12 MeV. It was concluded that 
a 10 MeV error is sufficient to cover the background uncertainty. The energies of the sim- 
ulated backgrounds were not varied during the fit. It was estimated that this introduces 
a negligible additional systematic uncertainty. 

3.5.5 Aspect ratio uncertainty 

As discussed in Section 2.5.3, uncertainties in the aspect ratio can bias the reconstruc- 
tion of jet angles and therefore the s' determination as well. 

Its effect on the centre-of-mass energy measurement was checked by introducing in 
simulation an artificial aspect ratio bias of a size of 3 x 10~ 4 . It consisted in scaling the 
longitudinal component of the momenta of charged particle tracks by this amount. The 
resulting change of the measured centre-of-mass energy was found to be 24 MeV. This 
value was taken as the aspect ratio uncertainty. 

3.5.6 Systematics from LEP beam parameters 

The effects of the Lep beam parameters on the energy measurement were taken to be 
the same as those found in the generator studies of dimuon events, namely an uncertainty 
of ±6 MeV on AE CM . 

3.5.7 Knowledge of the Z mass 

As in the dimuon analysis, there will be an uncertainty in the estimate of AE CM of 
±4 MeV, due to the uncertainty of the Z mass as measured [13] at LepI. 

3.5.8 Error summary 

Table 2 summarizes all contributions to the error of the measurement. The total 
systematic uncertainty on AE CM is 126 MeV. Thus the e + e~ — > qq(7) events gave a 
value of: 

AE CM = -0.116 ± 0.106 (stat.) ± 0.126 (syst.) GeV. 



4 Conclusions 

High-energy data collected in Delphi during the years 1997-2000 were analysed in 
order to cross-check the Lep energy measurement. The average difference between the 
centre-of-mass energy measured from radiative return e + e~ — > fi + fi~ (7) events and the en- 
ergy provided by the Lep Energy Working Group was found to be ±0.241±0.150 (stat.)± 
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QOUICe 


H/iioi on zau/cm ^iviev j 


Data/MC disagreement 


120 


ISR modelling 


p. 




T^vpi ornpnt p\ t ion 

1 X CllilllCllliClllWll 


2D 


Background 


10 


Aspect ratio 


24 


Lep parameters 


6 


MC statistics 


14 


Z mass 


4 


Total systematic error 


126 


Statistical error 


106 



Table 2: Contributions to the uncertainty on the centre-of-mass energy shift, as measured 
from radiative hadronic events. 

0.038 (syst.) GeV. The corresponding analysis using the e + e~ — > qq(7) events yielded a 
difference of — 0.116±0.106 (stat.) ±0.126 (syst.) GeV. Taking account of correlated sys- 
tematic errors (aspect ratio, Lep parameters, Z mass), the two results can be combined 
to give an overall Delphi result: 

AE CM = +0.073 ± 0.094 (stat.) ± 0.065 (syst.) GeV. 

Thus the Delphi data are compatible with the values reported by the Lep Energy 
Working Group [2], who quote values of the LEP centre-of-mass energy with uncertainties 
ranging from ±20 to ±40 MeV, depending on the year of operation. 
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